Draft version January 13, 2013 

Preprint typeset using I^Tf^X style emulateapj v. 08/22/09 



O 

m 



ON THE NATURE OF SPECTRAL LINE BROADENING IN SOLAR CORONAL DIMMINGS 

L. R. DOLLA 1 AND A. N. ZHUKOV 1 ' 2 
Draft version January 13, 2013 

ABSTRACT 

We analyze the profiles of iron emission lines observed in solar coronal dimmings associated with 
coronal mass ejections, using the EUV Imaging Spectrometer on board Hinode. We quantify line 
profile distortions with empirical coefficients (asymmetry and peakedness) that compare the fitted 
Gaussian to the data. We find that the apparent line broadenings reported in previous studies are 
likely to be caused by inhomogeneities of flow velocities along the line of sight, or at scales smaller 
than the resolution scale, or by velocity fluctuations during the exposure time. The increase in the 
amplitude of Alfven waves cannot, alone, explain the observed features. A double- Gaussian fit of the 
line profiles shows that, both for dimmings and active region loops, one component is nearly at rest 
while the second component presents a larger Doppler shift than that derived from a single- Gaussian 
fit. 

Subject headings: line: profiles — Sun: corona — Sun: coronal mass ejections (CMEs) 
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1. INTRODUCTION 

Coronal mass ejections (CMEs) are one of the most 
powerful energy release phenomena in the solar atmo- 
sphere. After such events, coronal dimmings usually ap- 
pear as regions of temporary lowered intensity in coronal 
lines, near the erupting active region. The fainter emis- 
sion is interpreted in terms of density decrease, as it is 
usually observed in several temperature regimes, and is 
usually linked to the coronal mas s removal that follows 
the eruption ([St erling & Hu dson! 119971 : lHarrison et al.l 
120031 : iZhukov fe Auch ere 2004 ). Spectroscopic studies 
showed that coronal dimmings are associated with large 
outflows with speed s of tens of kilometers per second 
lasting several hours ([Harra fc Sterl ing 200ll lHarra et all 
20071) ^ Using the EUV Imaging Spectrometer (EIS, 
Culhane et al. I [20071 ) on board the Hinode spacecraft, 



many observations revealed that strong outflows occur- 
ring in low-intensity regions were also associated with 
larger line width s than in the su r round ing area (e.g. 
iDel Zannal I2008T ). iDoschek et all (|2008f ) suggest that 
these wide profiles may be due to turbulence or to the 
addition of several components having different intrinsic 
velocities. They even note that some profiles show sec- 
ondary components that make the pr ofile depart from 
a p ure Gaussian profile. A ccording to iMcIntoshl (|2009[ ) 
and iMcIntosh et "ahl ([2009), the broadening is due to the 
growth of Alfven wave amplitude in the magnetically 
open and rarefied region of a dimming. All these stud- 
ies stress the possible connection between the solar wind 
and the observed outflows. 

From a theoretical point of view, both interpretations 
are possible. The corona is optically thin for most coro- 
nal emission lines. The observed lines are then the result 
of integration along the line of sight (LOS) and on scales 
smaller than the spatial resolution scales. This includes 
different layers having temperatures close to that of max- 
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imum formation of the line and also different magnetic 
flux tubes having differing plasma parameters. But in- 
tegration is also done over the exposure time. There- 
fore, Alfven waves and inhomogeneities of flows both 
produce spatial and temporal variations of Doppler ve- 
locities, which results in an apparent broadening of the 
integrated profile. The purpose of this paper is to in- 
vestigate the distortion of the line profiles to distinguish 
between both interpretations. 

2. ESTIMATING THE DISTORTION OF LINE PROFILES 

Different methods can be used to quantify the pres- 
ence of additional compon ents in line profiles. For ex- 
ample, llmada et al.l (|2008[ ) use a coefficient they call 
"additional component contribution to the line broad- 
ening", that makes use of the coefficients of a two- 
Gaussian fit that they apply in cases where the reduced 
chi-s quare (y 2 ) of the single-compon ent fit is larger than 
4.5. IMcIntosh fc De Pontieul (|2QQ9f ) evaluate the asym- 
metry by first finding the center of the line with a Gaus- 
sian fit, then calculating the difference of counts in sym- 
metric intervals in the red and blue wings, for a given 
offset velocity interval. 

The classical skewness and kurtosis of a statistical dis- 
tribution are influenced too much by discretization ef- 
fects in the line profiles we observe, while we prefer not 
to interpolate into the profiles to avoid adding informa- 
tion that would not be present before. They nevertheless 
inspired us to propose two empirical coefficients that we 
call hereafter "asymmetry" and "peakedness" to avoid 
confusion. They enable us to analyze the distortion of 
line profiles when the separation between the compo- 
nents is hardly noticeable, for cases where x 2 ~ 1 with 
a single-Gaussian fit. As for a x 2 computation, we use 
the squared differences between the real data and a fit- 
ted Gaussian. But contributions are positive or negative 
according to the sign of the difference and the data point 
position on predefined intervals of the line profile. 

Each coefficient C is defined as follows: 
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s k (X) ~ fk(X) 
ft (A) 



(1) 



where N is the total number of points where the con- 
tribution factor e(A), defined below, is non-zero. The 
spectrum Sfc(A) at wavelength A is discretized on spec- 
tral pixel fc, /fc(A) is the fit to Sfc(A), and Qe(A) is the error 
on Sfe(A). For the asymmetry, e(A) is defined as follows, 
with Ao being the center and a being the half-width at 
1 / y/e of the fitted Gaussian: 
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-1 if A G [A - 1.5 a; A -0.5 a] 
1 if A e (Ao - 0.5 a; A + 0.5 a) 

-1 if A G [A + 0.5a; A + 1.5 a] 
otherwise. 
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Tails beyond 2 and 1.5a, respectively, are not taken 
into account because they are more sensitive to noise. In 
practice, a coefficient equal to 1 means that over the main 
part of the profile the deviation was on average equal 
to one error bar, but systematically above or below the 
fitted Gaussian according to the above defined intervals. 
This can then hardly be a coincidence and appears as a 
significant distortion. 

The basic idea is that any departure from a Gaussian 
profile may be due to the presence of additional com- 
ponents having different centers, widths, or amplitudes 
at the line center. For simplicity, we further show syn- 
thetic profiles compound of only two components and 
fitted them with single Gaussians. But the results can 
be extrapolated to the addition of several components, 
with the limitation that any distortion can be smoothed 
out when their number increases. As a simple example, 
we simulated spectra with a dominant component having 
an amplitude at the line center of 100 in photon units, 
varying the relative amplitude of the other component. 
Artificial error bars are used for the consistency of the 
definition of the coefficients, using Poisson statistic (note 
that the coefficients become less sensitive as the statistics 
decrease, because of the increased error bars). But for 
the basic behaviors analyzed here, no real noise is added: 
it can create an artificial distortion, or on the contrary 
may blur the real one, and requires statistical analysis 
beyond the scope of this paper. 

First, we study the case of two components having 
identical widths a (taken equal to 1.4 pixel, the aver- 
age one measured in our data set, see below), with sep- 
arations of the centers varying from to 2 a; above 2 a, 
the double component becomes clearly visible. The co- 
efficient of asymmetry increases (in an absolute value) 
when the separation between both lines increases, but 
reaches its maximum when the ratio between their am- 
plitudes is around 50%. This is easily understandable: if 
it is too weak, the additional component cannot signifi- 
cantly distort the profile, and if both amplitudes are too 
similar, the profile is only broadened without showing 
any asymmetry. 



The upper panels of Fig. [T] show two examples for a 
second component of 50% amplitude, with a separation 
equal to a in panel (a) and 2 a in panel (b). In the case 
of panel (a), it is impossible to say by eye that there 
are two components, and the coefficients of asymmetry 
and peakedness (shown in the inset of each panel) are 
roughly equal to that is the value for a pure Gaussian 
(the first significant digit, not shown, is the third one 
after the decimal point). Nevertheless, the width of the 
resulting profile is already increased by 10% as compared 
to that of the components. In the case of panel (b), the 
asymmetry is difficult to notice without the visual com- 
parison with the fit. These two examples show that the 
fit tends to balance between the major component and 
the small tail by shifting the center of the fitted Gaussian 
from the major component toward the smaller one. The 
coefficients in Eqs. $Tj) to ([3]) were exactly designed to em- 
phasize this behavior. Note that while y 2 on ly estimates 
the deviation from the fitted Gaussian, the coefficient of 
asymmetry provides more refined information on the line 
distortion. Indeed, its sign gives the side on which the 
asymmetry is present: if positive, the profile has a major 
component on the shorter wavelengths (hence a tail at 
larger wavelengths), and vice versa. It can also detect 
asymmetries for the cases of x 2 values that would not be 
regarded as abnormal, i.e., close to 1. 

The coefficient of peakedness is less useful, because it 
is also influenced by the asymmetry of the profile. Never- 
theless, it is noticeably negative (i.e., the profile is flatter 
than a Gaussian) when the components are similar in 
intensity and separated by more than la (flat- top pro- 
file, Fig. [T] (c)). Therefore, it is sometimes possible to 
detect the presence of two (or more) components when 
the asymmetry fails to do so. The peakedness is pos- 
itive (i.e. the profile is more peaked) when the small 
component is less than 40% in amplitude than the major 
one, still for large separations. It is also positive when 
both components have different widths (Fig. [T] (d)). In 
this last case, as long as the components have the same 
center, the asymmetry stays around 0. 

3. OBSERVATIONS AND DATA PROCESSING 

We selected three Hinode /EIS data sets that show 
coronal dimmings associated with solar eruptive events 
(Table E]). We first concentrate on the Fe XII 195.12 A 
line because its high signal is suitable for the study of 
line profiles. 

3.1. Data set 1 : 2006 December 12-13 

In data set 1, NO A A AR 10930 was observed on 2006 
December 12-13 with exposure times of 30 s for every 
Y-position. This region presents bright structures in the 
first EIS raster (set la, see the first column in Fig. [2j) 
diverging roughly from X = — 120", Y = — 120". Some 
field lines are clearly connected to the core of the active 
region. There are no dat a from the Extreme ultravi olet 
Imaging Telescope (EIT. [Delaboudiniere et al1ll995[ ) on 
board the Solar and Heliospheric Observatory (SOHO, 
iDomingo etHI Il995[ ) before 23:49:26 on December 12, 
but our next EIS raster (lb, second column in Fig. [2]) 
shows that this area presents a dimmed region a couple 
of hours after the GOES B7.7 flare has started at 16:45 
UT. An X3.4 flare then occurred at 02:14 UT and was 
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Fig. 1. — Simulated profiles (asterisks) compound of two Gaussians (black solid lines) and fitted with a single Gaussian (gray solid line). 
(a): Sum of two Gaussians with the same width a = 1.4 pixels, and centers separated by la. The first one has an amplitude of 100 
counts, and the second of 50 counts, (b): same as a), but with a separation of 2a. (c): same as panel b), but the two Gaussians now have 
the same amplitude of 100 counts, (d): sum of two Gaussians of the same amplitude and same center, but widths of 1.4 and 2.8 pixels, 
respectively. The values in the insets are the coefficients of asymmetry and peakedness, respectively. The "0", "+" and "— " on top of the 
curves indicate the sign of the contribution e(A) in the intervals delimited by the vertical lines; dotted lines for the intervals associated 
with the coefficient of asymmetry (cf Eq. [2} and dashed lines for that of the peakedness (cf Eq. [3] note that we separately represent the 
intervals for the asymmetry and the peakedness in the first and the second row, respectively, for better visibility). Points that are more 
than 2a away from the line center are then not taken into account. 



TABLE 1 

EIS DATA SETS USED IN THIS STUDY. 



No. Date Start time a End time r b Comment 



la 2006 Dec 12 01:07:20 05:40:31 30 

b 12 19:07:20 23:45:45 30 After a B7.7 flare (16:45 c ) 

c 13 01:12:12 05:41:09 30 X3.4 flare at 02:14 

2a 2006 Dec 14 19:20:12 21:34:24 30 

b 15 01:15:19 03:29:31 30 After an XI. 5 flare (21:07 on December 14) 

c 15 04:10:12 06:24:24 30 Dimming partial recovery 

3a 2007 Aug 22 13:33:46 17:55:56 60 After a B1.2 flare (11:29) 

b 23 01:55:43 06:17:53 60 Overlap with an eruption 



a All times are UT. 

b Exposure time is in seconds. 

c All the flare times are flare start times as measured by GOES. 
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also associated to the dimming in the same area. This 
event was observed by EIT. The EIS slit started to raster 
this dimming area one hour after the flare b egan (set lc, 
see Fig. [3l also analyzed bv lAsai et al.ll2008l ). All rasters 
are scanned from west to east. For the Doppler refer- 
ence wavelength in raster lc, we only used the southern 
"quiet" part (pixel rows 15 to 40) to reduce biases due 
to the large Doppler shifts present on most of the map. 

3.2. Data set 2 : 2006 December 14-15 

Data set 2 corresponds to the same NOAA AR 10930 
observed several hours after data set 1, again with 30 s 
exposure times. Set 2a (first column in Fig. |4]) shows 
that the dimmings have recovered since the X3.4 erupt- 
ing flare. An XI. 5 fla re started at 21:07 UT on Decem- 
ber 14 . As noted by lHarra et all (|2QQ7l ) and iMcIntoshl 
(j2009[ ). a new dimming appears at the same location as 
during the December 13 event (second co lumn, set 2b). 
Such a kind of repetition was described bv lChertok et al.l 
(|2004f ). After a few hours, the dimming started to recover 
(set 2c). 

3.3. Data set 3 : 2007 August 22-23 

Data set 3 (Fig. [5]) corresponds to an active region 
observed on 2 07 Au gust 22 and 23, also analyzed by 
iDoschek et al.l ([2008) . For these two rasters, the expo- 
sure times were 60 s. There is a B1.2 flare 2 hr before 
the raster 3a, associated with dimmings that are best 
seen on images taken by the Extreme Ultraviolet Im- 
ager (EUVI), which is a part of the Sun Earth Connec- 
tion Coronal and Heliosp heric Investigation (SECCHI) 
suite (jHoward et a l. 2008) on board the Sola r TErrestrial 
REla tions Observatory (STEREO) mission (jKaiser et al.l 
2008). The difference images in Fig. [6] show a dimming 
in the southeast part of the active region, and a bright 
fanning-out structure that is preserved all along the ob- 
servations. 

A second event occurs just before the beginning of the 
second raster (data set 3b). A faint dimming appears 
and the loop system clearly produces an eruption from 
01:46 UT to 02:46 UT (Fig. [7j), spanning the fanning- 
out structure. Around 03:56 UT (last panel), the scan is 
in the middle of the raster and the dimming area does 
not increase any more (with even some faint recovery), 
while bright post-eruptive loops have started to be seen 
in the active region core. There is no recorded flare in 
the GOES catalog corresponding to this event, although 
one can observe a brightening in the GOES light curve 
around 3:00, that corresponds to the appearance of the 
bright loops in the EUVI images. 

All three data sets exhibit the same properties that 
will be discussed in Sect. HJ 

3.4. Data processing 

By fitting the line profiles with single Gaussians in- 
cluding a constant and a linear term, we obtain the peak 
intensity of the line, the LOS Doppler velocity and the 
non- t hermal velocity deri ved from the line width (see, 
e.g., iDere fc Mason| [T993). For the Doppler shift, abso- 
lute wavelength calibration is difficult, so like in many 
other studies, we take an average wavelength position 
along the Y-direction as a reference, after correcting for 
the slit tilt and orbital variation (standard procedures in- 



cluded in SolarSoft). Although we will show that the con- 
cept of non-thermal velocity has to be taken with care, we 
derive it by assuming an ion temperature equal to that of 
maximum line formation. An in strumental width equ iv- 
alent to 2.5 pixels is subtracted ([Doschek et al.ll2008f ). 

When applying the eis_prep.pro SolarSoft proce- 
dure, we used the recommended option for interpolation 
of the missing data (e.g., due to the cosmic-ray removal). 
This procedure leaves missing data empty when interpo- 
lation has no sense (more than three neighboring pixels in 
the Y-direction missing). If such missing data remained 
in the central half of the spectral window, the derived 
Doppler velocity, width or coefficient of distortion would 
be treated as missing data and set to in our plots. For 
aesthetic reasons, the corresponding pixels in the inten- 
sity images are replaced by a median of the neighboring 
pixels. 

We use the pointing coordinates given in the EIS fits 
files: no absolute pointing precision is necessary for this 
study. It is well known that the poi nting varies with 
wavelength (jYoung et al.ll2007bL I2009I ) 3 . When making 
comparison between different lines, we realigned the im- 
ages using the Fe XII as a reference. The misalignment 
in the Y-direction is corrected by using the procedure 
eis_ccd_of f set .pro in SolarSoft. There is a 2-pixel 
difference in the X-direction between images taken with 
the short wavelength detector (170 — 211 A) and the long 
wavelength one (246 - 292 A). 

4. RESULTS 

As shown in Fig. [2] to [5j bright loops in the core of the 
post-flare acti ve region present large redshifts in bot h 
legs (see also iWinebarger et all 120021 : iDel Zannal [20081 ). 
while the dimming regions present large blueshifts that 
can exceed 50 kins" 1 in what appears to be open mag- 
netic field lines or very lo ng loops. Both kinds of regions 
show large line widths (see Doschek et al. 2008). In set lc 
(Fig. [3]), the dimming area that presents large blueshifts 
and large widths extends out of the low-intensity part, 
into the quiet-Sun loops region. This area is nevertheless 
consistent with the area covere d by the dimming in the 
EIT difference images (see, e.g. lAsai et all 120081 , their 
Fig. 6). 

4.1. Asymmetric profiles of the Fe XII 195 A line 

Regarding the asymmetry and peakedness of the Fe XII 
195 A line, we note large areas of distorted line pro- 
files, appearing in black or white in the last two rows 
in Fig. [2] to [5] The most noticeable result is that neg- 
ative or positive asymmetries cluster and do not spread 
around randomly. The largest distortions appear in areas 
where the non-thermal velocity is very large. Asymme- 
tries are predominantly positive in bright areas (cooling 
post-flare loops), while negative coefficients of asymme- 
try are essentially present in the low-intensity areas (dim- 
mings). Globally, red(blue)shifted areas correspond to 
bright (dark) areas. 

The patterns in peakedness are less clear. Globally, 
bright loop areas present a positive peakedness, i.e., pro- 
files more peaked than a pure Gaussian. There are some 

3 see also http : //msslxr . mssl . ucl .ac.uk: 8080/eiswiki/Wiki . j sp?page=CCDC 
and http : //msslxr . mssl . ucl .ac.uk: 8080/eiswiki/Wiki . j sp?page=CCD0f f setX 
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Fig. 2. — From top to bottom: Fe XII 195.12 A line peak intensity, Doppler velocity, and non-thermal velocity observed by Hinode /EIS 
on 2006 December 12 and 13 (data sets la and lb). Also shown are the derived coefficients of asymmetry and peakedness. 
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Fig. 3. — Same as Fig. [2] but for data set lc. Red, black and 
blue contours, in their respective panels, show the areas with a 
non-thermal velocity larger than 60 kms -1 . 



exceptions, like the negative areas in the core of the ac- 
tive region in set lc (Fig. [3]). But a closer look at the 
spectral profiles showed us that they correspond to large 
asymmetries, so that the fit is heavily shifted from an 
otherwise not particularly flat profile; this weights the 
coefficient of peakedness and shows that this coefficient 
can be better defined in the future. 

Some examples of spectra extracted from the dimming 
area in data set lc that presents large non- thermal ve- 
locity (contoured area in Fig. [3]) are shown in Fig. [8j for 
different combinations of both coefficients. Note that in 
most cases, the reduced x 2 is relatively small, especially 
when statistics are poor, even though the distortion of 
the line profile can be detected by the eye. Our coeffi- 
cients better detect the mismatch between the line profile 
and the fitted Gaussian. Following our analysis in Sect.[2j 
we interpret these asymmetries as due to the presence of 
a minor component in the profiles: in the blue (red) wing 
for the negative (positive) asymmetry. This is consistent 
with multiple flows merging in the same LOS or pixel. 

We checked the results produced when not using the 
interpolation procedure for the missing data, and the 
asymmetry and peakedness show the same overall behav- 
ior. Therefore, the effect of the interpolation procedure 
fo r explaining the o bserved features can be ruled out. 

lHara et all ([20081 ) also found, in upflows, deviations 
from a single Gaussian in the blue wing of the Fe XIV 
274 A and the Fe XV 284 A lines. Note that we con- 
sider here smaller outflows and mu ch less obviou s asym - 
metries than those pr esented by llmada et al.l ([20071 ) , 
llmada et all ([2008D or lAsai et al.l (j2008|) for the same 
data set (lc), even if we also detect the prominent fea- 
tures that they found. In particular, the small black spot 
near X = 360 ", Y -90 " in F ig. [3] corresponds to the 
BS1 region of I Asai et all ([2008D . The method we pro- 
pose enables us to emphasize more subtle distortions, and 
shows that they are widely present in solar active regions, 
including parts with downflows. In the contoured area 
of data set lc, i.e., in pixels with non-thermal velocities 
larger than 60 kms -1 (Fig. [3]), we found 21% of pixels 
having an absolute coefficient of asymmetry larger than 
0.5. This goes up to 28% if we count pixels with asymme- 
try or peakedness larger than 0.5. In bright loops, most 
of the pixels present noticeable distortion of the spectral 
profiles. 

4.2. Double- Gaussian fit 

For the multi-component analysis, we make double- 
Gaussian fits to the line profiles (see details in Ap- 
pendix [A]). Both components of the double- Gaussian fit 
for data set lc are over-plotted in the profiles shown in 
Fig. [H In Fig. [9] we present Doppler shifts associated 
with both components of the Fe XII 195 A line (top row) 
for data se t lc . Besides pixels with missing spectral data 
(see Sect. 13.41 and e.g., black pixels in the third panel 
of Fig. [3]), we rejected pixels where the components are 
separated by more than 2 Gaussian widths of the first 
component. This criterion is used to prevent us from 
considering cases where the double fit selected neighbor- 
ing lines of different rest wavelengths, that appear in the 
195.12 A spectral window in some pixels. We also re- 
jected pixels where the components are separated by less 
than 10 kms -1 (nearly one third of a spectral pixel), or 
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Fig. 4. — Same as Fig. [2]but for data sets 2a, 2b, 2c, observed by Hinode/ElS on 2006 December 14 and 15. 



8 



L. Dolla & A. Zhukov 



2007/08/22 13:33:46 

Peak intensity 



2007/08/23 01:55:43 

Peak intensity 




-700 -650 -600 -550 



Doppler Velocity 




-700 -650 -600 -550 



Non-thermal velocity 




-700 -650 -600 -550 



Asymmetry 




-700 -650 -600 -550 



Peakedness 



-150 

-200 
-250 
-300 



Ml m 




-600 -550 -500 -450 



Doppler Velocity 




-600 -550 -500 -450 



Non-thermal velocity 




-600 -550 -500 -450 



Asymmetry 



-100 
-150 

-200 




-600 -550 -500 -450 



Peakedness 



-100 




, i , , , , i , , , , 




-150 




iw 


'. 


-200 






I 




-250 










-300 











-700 -650 -600 -550 

Solar X 



-600 -550 -500 -450 

(arcsec) 



1.0 
0.5 



-0.5 
-1.0 



1.0 

0.5 

0.0 
1-0.5 
1-1.0 



Fig. 5. — Same as Fig. [2] but for data sets 3a and 3b, observed 
by Hinode /EIS on 2007 August 22 and 23. 



where the peak intensity of one of the components is less 
than 20% of the other. We applied these conservative 
^filters to ensure that the double component that is fur- 
ther analyzed is significant. The white pixels in the right 
\ column of Fig. [9] mostly correspond to the rejected pix- 
£ els. Judging by the retained pixels in this figure, most of 
P the faster components are well separated from the slower 
^ones, with a contribution of more than 20% for the mi- 
gmor component, both in the dimming area and the active 
region core. Because of the limitations of the double- 
Gaussian fit explained in Appendix [Aj a thorough anal- 
ysis of the ratios of amplitude between the components 
is left for future work. 

The first row of Fig. [10] shows the distribution of 
Doppler velocities of both components compared to that 
of the single- Gaussian fit, for pixels where the non- 
M thermal velocity is larger than 60 kins" 1 (i.e., the con- 
| toured areas in Fig. [3]). We plot the distributions sepa- 
rately for pixels in the dimming area and for the active 
region core. In the active region core, we only retained 
pixels where the single- Gaussian fit provided a redshifted 
component. This prevents us from including the highly 
blueshifted event related to the flare that is studied by 
lAsai et aTl ([20 08), and focuses the analysis on the cool- 
ing loops. The medians of these distributions are listed 
in Table [2] As the error on the velocity determination is 
around a few kilometers per second, we can conclude that 
the slower component, in the dimming area and the ac- 
tive region core, corresponds to emission from a roughly 
static plasma. For this reason, we will hereafter call this 
component the (nearly) static component. The second 
component presents Doppler velocities larger than those 
derived from the single- Gaussian fit. We will call it the 
dynamic component. 

Regarding the non-thermal velocities retrieved from 
the component widths in the dimming, we find distri- 
butions with medians at 70, 42, and 76 kms -1 for the 
single- Gaussian fit, the static component, and the dy- 
namic component, respectively. Therefore, the double 
component analysis shows that the static component 
presents a non-thermal velocity more akin to what can 
be found in the quiet Sun, even though it may consist 
of more components, so that the individual profiles may 
get even narrower. 

The reliability of the double-component fit was tested 
on synthetic line profiles (see Appendix [A]). This analysis 
shows that: (1) the presence of two components in the 
data is real and is not an artifact of simply applying 
the double-component model on the broad profiles; (2) 
the average non-zero value of the static component is an 
artifact of the fitting procedure and cannot be considered 
as a definite result. 

Overall, the analysis we made in Sect. [2] is verified: 
areas with a definite positive (negative) asymmetry cor- 
respond to a redshifted (blueshifted) dynamic component 
in the profile. We also note that the redshifted dynamic 
component in the loop areas is larger in the core of the 
active region than in its surroundings. 

4.3. Results for other spectral lines 

Several other lines were recorded in our data sets. We 
concentrate again on data set lc that shows the asym- 
metric profiles more clearly. As the signal in other lines 
is lower, we spatially rebinned the data by 3 pixels in the 
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Fig. 6. — Leftmost 195 A STEREO-SECCHI/EUVI panel (spacecraft B) shows the last pre-emption image of the active region observed 
in data set 3. The other panels are base difference images with the image in the first panel subtracted from other images, taking into 
account the differential rotation. The images in the third and fourth panels were taken at the start and at the end of the raster of data 
set 3a. 
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Fig. 7. — The leftmost 195 A STEREO-SECCHI/EUVI panel (spacecraft B) shows the active region observed in data set 3, 50 minutes 
before the beginning of raster 3b. The other panels are base difference images with the image in the first panel subtracted from other 
images, taking into account the differential rotation. In the third panel, one can see dimmings and a loop system erupting. In the last 
panel, which corresponds to the middle of raster 3b, one can see bright loops appearing. 



TABLE 2 

Medians of the distributions of velocities for data set lc a . 







Dimming 


area 


Active region 


core b 


line ID 


^single 


^static 


^dynamic 


^single ^static 


^dynamic 


Fe XII 195.12 A 


-27 


-6 


-50 


24 1 


64 


Fe XIII 202.04 A 


-29 


-4 


-72 


10 -2 


30 



a corresponding to the contoured areas in Fig. [9] i.e., non-thermal velocities larger 

than 60 kms -1 with the single-Gaussian fit. All velocities are in kms -1 . 

b For the active region core, only pixels with redshifted single-Gaussian component 

are retained (see the text). 

c Derived from the single-Gaussian fit. 

d Derived from the double-Gaussian fit. 



X-direction and 3 pixels in the Y-direction. We verified 
that this rebinning does not produce additional distor- 
tions (Appendix [B]). This also demonstrates that the 
source of the distortion has nothing to do with the res- 
olution scale: it has to be found in the LOS integration, 
in scales much smaller than the resolution scale or in 
temporal variation. 

After the rebinning, we derived the coefficient of asym- 
metry and peakedness for several lines: Fe XII 195.12 A, 
Fe XIII 202.04 A and Fe XIV 274.20 A (see Fig. [TT]) . 
These ions have peak formation temperature increasing 
from 1.6 to 2 MK. As there is a wavelength-dependent 
shift of the spatial position on the detector, we co-aligned 
all maps with that of the Fe XII line (see Sect. 13. 4|) . To 



ease the comparison, we overplot the 60 kms -1 contours 
for the non-thermal velocity of Fe XII ions on all the pan- 
els. Different thresholds for the color scale are used to 
ease the visual comparison. It must be emphasized that 
the differences in the signal statistics greatly influence 
the value of the coefficients due to the inclusion of the 
error bars in Eq. [IJ 

Areas of large line widths (blue contours) appear more 
fragmented when the formation temperature increases, 
but this may be a selection effect due to the definition of 
the non-thermal velocity as the excess width as compared 
to the formation temperature of the line. The contours 
are mainly co-spatial for different lines in the dimming 
area. We also note that the same areas of negative or 
positive asymmetry or peakedness appear at the same 
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Fig. 8. — Sample of Fe XII 195.12 A line profiles in data set lc. We selected pixels in the dimming area (X< 250") where non-thermal 
velocities are larger than 60 kms - 1 (contoured area in Fig.[3j), for combinations of distortion coefficients below —0.5, between —0.5 and 0.5, 
or above 0.5. The inset in each panel contains the coefficients of asymmetry and peakedness, in this order. The thick gray line represents the 
best single-Gaussian fit to the observed profile, with the associated x 2 mentioned in each panel. Vertical lines indicate the rest wavelength. 
The thin gray lines represent the two components of the double-Gaussian fit. 



locations for every wavelength, especially in the dimming 
area for Fe XII and Fe XIII. We do not show the Doppler 
and non-thermal velocity maps, but they are similar to 
those taken in the Fe XII line, as has already been show n 
in several studies (e.g. Im ada et al.ll2007l : lJin et al.ll2QQ9h . 
No large distortions are found in the dark pixels of the 
central dimming area for Fe XIV, which is essentially due 
to the lack of photon statistics there. 

Among the other available lines, the Fe VIII 185.21 A 
line is too faint and leads to very noisy maps of asym- 
metry, peakedness and even line widths. For the Fe X 
184.54 A line, we find distorted profiles in the active re- 
gion, but the coefficients of distortion are close to in the 
dimming area, even though we could discern some black 
patches reminiscent of what we find in other wavelengths. 
Again, this may be a question of lack of photon statis- 
tics. The Fe XI 188.23 A line is composed of two close 
lines of the same ion. This would make the distortion 
analysis difficult. The analysis of the Fe XV 284.16 A 
line is complicated by blending issues. 

We also find a similar behavior for data set 3. Positive 
(negative) asymmetry is seen in the loop (fanning-out 
structure) area, for the same lines as shown in Fig. [11] 
and for additional lines, including another Fe XII line 
at 193.51 A and another Fe XIV line at 264.78 A (not 
shown) . 



As demonstrated in Appendix [Aj the double-Gaussian 
fit of closely situated lines is not precise enough to allow 
a reliable comparison pixel by pixel. This is especially 
true for different wavelengths, for which one has also to 
correct accurately pixel shifts on the detector due to the 
wavelength dependence. This effect is even more impor- 
tant when the lines are recorded on the different EIS 
detectors. We show the Doppler maps and velocity dis- 
tributions for the Fe XIII 202 A line in Fig. [9] and Fig.flOl 
respectively. 

5. DISCUSSION 

We found a clear pattern of asymmetry of coronal line 
profiles. In dimmings, a substantial part of them is 
left-skewed (negative asymmetry), while in bright loops 
most of them are right-skewed (positive asymmetry). 
The largest distortions appear in areas where the non- 
thermal velocity or, more precisely, the width of the 
profile, is very large. This is especially true in the 
left-skewed profiles. A double- Gaussian fit reveals that 
these asymmetric profiles are composed of one nearly 
static component and (at least) one dynamic component. 
Most of the time the dynamic component is systemati- 
cally blueshifted (redshifted) in the case of a left-skewed 
(right-skewed) profile. 

5.1. Alfven waves or flow inhomogeneities ? 
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Fig. 9. — Slower component (or "nearly static", left column) and faster component (or "dynamic", right column) of the double-Gaussian 
fit made for data set lc. The top row corresponds to the Fe XII 195 A line and the bottom row to the Fe XIII 202 A line. White pixels in 
the faster component panels mostly correspond to fits not taken into account for analysis (see the text for details). 

duration events), may be another explanation. It is note- 
worthy that these interpretations are compatible with the 
fact that a lot of line profiles do not show any substantial 
distortion in the dimming areas that nevertheless present 
a broad profile (e.g., the central panel of Fig.[8j). Indeed, 
for small differences in the centroid positions of the com- 
ponents or for small relative intensity of the additional 
component, the resulting profile is broadened but does 
not present a noticeable distortion, like in Fig. [T] (a). In 
this simple two-component case, the width of the result- 
ing profile is increased by about 10% as compared to that 
of the single components. This turns into an increase of 
30% in non-thermal velocity in the case of an Fe XII line. 
This value is larger than the average increase observed 
in the dimming area of data set 2 as shown in Fig. [4] (see 
also p anel (D) in the Figure 3 of the work by IMcIntoshl 
120091 ) . The lack of photon statistics in dimmings may be 
another reason for non-detected distortions. But the ad- 
ditional component is probably present in all pixels with 
large line profiles, and thus can explain the observed in- 
crease of the coronal line widths. 



It appears that transverse velocity components of 
M HD waves and, in particular, Alfven waves as suggested 
by IMcIntoshl (|2009l ) cannot be responsible for the line 
broadening observed in dimmings. If this were true, the 
broadened profiles should be symmetric, because veloc- 
ity perturbation in a wave oscillates around zero. If one 
assumes that the spectrometer takes a snapshot of MHD 
waves having periods longer than the exposure time (30 s 
per X-position), one should find an approximately equal 
number of points with negative or positive asymmetry, 
randomly positioned. One could also imagine that all 
the waves in adjacent pixels are in phase so that their 
velocity perturbations all have the same sign. But this 
would have to be the case not only in the Y-direction (si- 
multaneous exposure) but also in the X-direction, which 
implies the coincidence of the rastering speed with the 
wave propagation speed. This appears very unlikely, es- 
pecially in so many different data sets. Above all, waves 
cannot explain the fact that the left-skewed profiles are 
found in dimmings and the right-skewed profiles in loops. 

Our interpretation appears more plausible: the broad- 
ened profiles are due to the additional component that 
has the same velocity direction as that found with the 
single- Gaussian fit. Downflows in loop s may correspon d 
to draining of cooling material (e.g. Bradshawl l2008f ) . 
while upflows correspond to material escaping the corona 
or refilling it after the eruption. The multi-component 
approach offers a new point of view on the "corona l cir- 
culation" concept suggested bv lMarsch et al.l (|2008f ): the 
dynamic component is most of the time observed to- 
gether with a nearly static one in the same LOS. Plasma 
velocity varying quickly on temporal scales smaller than 
the exposure time (associated with jets or any other short 



5.2. Effect of line blending 

We now assess the effect of line blending on the ob- 
served distortions. First of all, let us remind that to 
calculate the coefficients of asymmetry we select only a 
window around 2a from the line center (see Eq. [2]). This 
means that one can consider any line lying beyond that 
interval as not contributing to this coefficient. The inter- 
val is even smaller for the peakedness (see Eq. |3j). Note 
also that when considering the double- Gaussian analysis, 
we discard any fits where the center of the second compo- 
nent is beyond 2a from the main component. This filter 
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Fig. 10. — Distribution of Doppler velocities measured in the pixels of data set lc where the non-thermal velocity in the Fe XII 195 A 
line is larger than 60 kms -1 (calculated using a single-Gaussian fit). Left column: for pixels in the dimming area (X < 250"). Right 
column: for the redshifted active region core (see the text). The top row corresponds to the Fe XII 195 A line and the bottom row to the 
Fe XIII 202 A line. The black histogram corresponds to the result of the single-Gaussian fit, the red one to the static component of the 
double-Gaussian fit (around 0) and the green one to the dynamic component. The bin size is 5 kms -1 . 

is applied after all the iterations of the fitting procedure 
(not during the iterations), which means that we have 
no selection effect. Note also that a very few of such fits 
were discarded. They are located mainly in the active 
region near X=300", Y=-120" in the case of Fe XII (see 
Fig. [9j where they appear as the white pixels that do not 
correspond to the black pixels of the line- width panel in 
Pig.©. 

A 2cr interval for the Fe XII 195.12 A line makes on 
average 0.06 A (or about 92 kms -1 ). It can go up to 
0.1 A in the large widths areas (dimming and active re- 
gion core). iBrown et al.l (|2008f ) report on the following 
neighboring lines for Fe XII 195.12 A. There is an Fe XII 
194.92 A line in the red wing. This is more than 4<j apart 
(and it is only less than 2% of the 195.12 A line, accord- 
ing to these authors). On the blue wing, the closer one is 



Fe XIV 195.25 A, which is more than 4<j apart from the 
195.12 A line, or 2.5a if we consider the large width case 
(and again, it is less than 2% of the intensity of the Fe XII 
line). Thus, this line is not taken into account in our 
computations and analy sis. In the CHIANTI database 
(jDere et al.l 119971 120091 ). we also find a Ni XVI line at 
195.27 A which is again too far. The only critical blend- 
ing issue is another Fe XII line at 195.18 A. For a den- 
sity of 10 10 cm -3 , i.e in active regi ons, it is predicted to 
contr ibute around 10% in intensity ([Del Zanna fe Masonl 
120051 ) . When we simulated such kind of line profiles with 
the high photon statistics characteristic of active regions, 
we found coefficients of asymmetry larger than 1. But 
the coefficients really measured in data set lc are much 
larger. Moreover, the distortion is hardly distinguishable 
by the eye in the simulated profiles, contrary to that ob- 
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Fig. 11. — Coefficients of asymmetry (left column) and peakedness (right column) for the Fe XII 195.12 A, Fe XIII 202.04 A and Fe XIV 
274.20 A lines. EIS pixels are binned by 3 in each direction. Red contours correspond to 60 kms -1 for the non-thermal velocity of Fe XII, 
while the blue contours correspond to the same limit for each of the two other lines. 



served in the real active region profiles. This suggests 
that the main contribution for the positive asymmetry 
of the Fe XII 195.12 A line profile in loop areas is not 
the self-blend with Fe XII 195.18 A. The self-blend can- 
not account anyway for the negative asymmetry in the 
dimming. 

According to lYoung et all (|2007a[ ). the Fe XIII 
202.04 A line is unblended, the Fe XIV 274.20 A line 
blend with a Si VII line (in the blue wing ) can be ne- 
glected in active regions (Lian get al.ll2010f ). 

Therefore, the negative asymmetry in dimmings, free 
of any blending, can be ascertained from the Fe XII 
195.12 A and Fe XIII 202.04 A lines. The positive asym- 
metry in loops, at least in the active region core, can 
be ascertained from the Fe XIII 202.04 A. This posi- 
tive asymmetry can also be ascertained from the Fe XIV 
274.20 A line because the weak blend it is subject to is 
present in the blue wing, which means it should produce 
a negative, and not positive asymmetry. Overall, simi- 
lar patterns of asymmetry and/or peakedness observed 



in different lines (see Fig. [TT]) suggest that the line pro- 
file distortions are not essentially produced by blending. 
Inhomogeneities of flow then appear as the main cause 
of these asymmetries, even though line blending compli- 
cates the analysis of the components. 

5.3. Flow geometry and link with eruptions 

We showed that rebinning the data does not produce 
additional distortion of line profiles. On the contrary, it 
helps to emphasize the pre-existing distortions. Conse- 
quently, the inhomogeneities in velocities on scales larger 
than the resolution scale are not sufficient to produce the 
observed distortions: such inhomogeneities must be al- 
ready present before the binning. They should exist on 
scales smaller than the initial resolution scale (~ 1"), or 
must be due to LOS or temporal integration. 

An interesting question is the nature of the static com- 
ponent, most of the time the more intense one. It could 
correspond to a zero sum between blue and red shifts 
that will symmetrically broaden the static component in 
excess of the thermal width (the usual interpretation of 
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the non-thermal velocity). This can be caused by waves 
or turbulence, for example Alfven waves. Another possi- 
ble interpretation is the presence of flux tubes where the 
plasma has no or small overall projected velocity on the 
LOS. In addition to a truly static component, this may 
also be the case when the LOS (or pixel) superposition of 
many flux tubes with different (counter-streaming) flow 
velocities is present. Additionally, the cross section of 
some of the flux tubes may be perpendicular to the LOS 
(e.g., near the loop apexes). These interpretations can 
hold both for loops (multiple parallel strands in one loop 
or overlying different loops) or the dimming areas. In 
the latter case, the areas of negative asymmetry seem to 
present a radial pattern (especially for data sets lb, lc, 
3a and 3b, see the black areas in the corresponding pan- 
els in Fig.[2j[3]and[5j). It is tempting to link them to flux 
tubes of the fanning-out structure. The static compo- 
nent could then correspond, here again, to neighboring 
radial strands with plasma at rest or to apexes of very 
long loops overlying the dimming area. In three dimen- 
sions, one can imagine a mix of field lines, some of them 
nearly vertical and maybe open, some other ones more 
inclined or even horizontal, all more or less entangled, 
with proportions depending on the local topology of the 
magnetic field. 

Note also that more components may be present in 
the line profiles. Due to their closeness, it is very dif- 
ficult to distinguish them. The profiles do not present 
enough "bumps" to unambiguously retrieve every com- 
ponent by a multi-component fitting procedure. With 
the present analysis, the only thing we can be certain 
of is that these profiles are not symmetric enough to be 
explained by a Gaussian or any other symmetric distri- 
bution of velocities for the components. There are at 
least two components. 

It is tempting to relate the systematic blueshift of the 
nearly static component to a real flow corresponding, 
e.g., to the nascent solar wind. But first, one has to 
remember that this Doppler velocity has been derived 
by taking a reference wavelength computed by using the 
centers of the single-Gaussian fits, i.e., not taking into ac- 
count the fact that some of these profiles have multiple 
components. Second, this systematic shift may be a bias 
of the two-Gaussian fit (see Appendix lAj) . Third, these 
Doppler shifts are within the error bars. This prevents 
us from firmly concluding on the reality of this upflow. 

Another interesting issue is the relation between neg- 
ative asymmetries and eruptions. In our data sets, it 
is obvious that the negative asymmetries get larger and 
more extended after every eruption (cf Sect. [3] and Fig. [2] 
to [5]). We note that data sets 1 and 2 are taken in the 
same active region that underwent a significant rotation 
from near the disk center toward the west limb. Never- 
theless, the Doppler shift patterns as shown in Fig. [2j [3] 
and H] look very similar. We could not find any marked 
difference between the flows due to the different positions 
on the disk. 

Three effects can explain the increase of asymmetries 
after eruptions. First, we can suppose that the high- 
speed dynamic component is already present before the 
eruption. The disappearance of a large part of the coro- 
nal material, especially in the apexes of closed loops dis- 
cussed above, may lower the intensity of the static com- 
ponent and thus make the dynamic one more prominent 



in the profile. This can be the case for data sets 1 and 
2. This interpretation is reinforced by the fact that be- 
tween data sets lc and 2a, the transient coronal hole has 
been refilled, while the negative asymmetries have dis- 
appeared. The case of data set 3 is more complex: the 
observed dimmings here are very faint and do not really 
correspond to the most blueshifted structure. 

Still with a pre-existing dynamic component, a second 
explanation is that the reorganization of the magnetic 
topology during the eruption leads to a better alignment 
of the flows with the LOS (e.g., more vertical flux tubes, 
as expected during a dimming event). This increases 
the LOS component of the flow and the asymmetry of 
the resulting profiles. This explanation can account for 
the transition between data sets 3a and 3b, or maybe 
for the transition from data set 2a to 2b. But there 
is a large area of redshifts in data set lb around X = 
50", Y = — 80". After the eruption, these areas show 
only blue shifts. This implies that the flows reverted from 
downward to upward. It is then difficult to envisage that 
all the outflows were present before. 

The third possibility is that the dynamic component 
may be a feature appearing in dimming only after erup- 
tions. It is possible that dimmings are magnetically open 
struc t ures (see e.g. iKahler fc Hudson! 120011 : lAttrill et al.l 
120081 : lMcIntosh|[2009l ). The plasma then can easily es- 
cape from the Sun, pr oducing the observed outflows (e.g. 
iMcIntosh et al.ll2010[ ). From this point of view, it is pos- 
sible that apparently broader profiles (cf Sect. 14. 2[) of the 
dynamic component correspond to the acceleration that 
produces sub-components of increasing blue shift during 
the exposure time. They would then be considered as 
a single dynamic component by the double- Gaussian fit. 
Therefore, the increase of asymmetry of the coronal line 
profiles may be due to an increase of the transient activ- 
ity after the eruption. Whatever the case may be, the 
outflow corresponding to the dynamic component may 
take part in the refilling of the corona during the dim- 
ming recovery. 

It is worthwhile to emphasize that most of the LOS 
velocities that we find for the dynamic component are 
subsonic (the sound speed c s w 200 kins" 1 at 1.6 xlO 6 K; 
cf Fig. [TO]). 

6. CONCLUSIONS 

We analyzed the profiles of coronal iron lines observed 
with Hinode/ElS by using empirical coefficients (asym- 
metry and peakedness) efficient at both detecting and 
discriminating different distortions of the line profiles. 
These coefficients are sensitive to the photon statistics in 
the profiles, but this is what makes them reliable. They 
are a good complement to multi-component analysis, 
particularly for a preliminary check of interesting data 
sets, as they are fast to compute. They are especially 
suitable to reveal the presence of very close components, 
which are difficult to resolve with a double-Gaussian fit. 

We showed that the line broadenings correlated with 
large Doppler shifts observed in several active regions, 
are associated with the presence of more than one Gaus- 
sian component in the line profile. Upflows are found in 
the coronal dimming area and downflows in the loops , 
confirming previous st udies (e.g. lHarra fc Sterling! 120011 : 
IWinebarger et a3l2002f ). However, their velocities are un- 
derestimated if a single- Gaussian fit for the line profile is 
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used. Both static plasma and large flows (several tens of 
kilometers per second) are present in the same LOS or 
pixel. 

We do not rule out the interpretation that the inherent 
density decrease in coronal dimmings can broaden the 
line profiles due to the induced increase of the A lfven 
wave amplitude, as sugg ested bv iMdntoshl (j2009). but 
this cannot be the sole cause. It is difficult to explain 
the asymmetries observed on large scales only by Alfven 
waves. One has to take into account the effect of inho- 
mogeneities of velocities in the LOS (overlapping strands 
with different inclination or intrinsic velocity) or in the 
spatial resolution pixel (flows with sub-resolution struc- 
ture) , or transient variations of the velocity on timescales 
smaller than the exposure time (30 s). The concept of 
non-thermal velocity derived from a single- Gaussian fit 
must then be handled with care. In particular for on-disk 
observations, one cannot directly interpret this quantity 
as the measure of average fluctuations of velocities in the 
observed pixel, especially in terms of Alfven wave am- 
plitude. We demonstrated that this pattern is visible in 
three types of coronal structures which exhibit large-scale 
flows: CME-associated dimmings, loops, and fanning-out 
structures. A new field of investigation is now opened to 
analyze the structure of inhomogeneities in an optically 
thin coronal medium. 

Finally, we note that since the submission of this paper, 
three more works emphasize d the importance of a sym- 
metric coronal line profiles. Mcintos h et al.l (j2010h use 
a different approach bas e d on the method described by 
iMcIntosh fc De Pontieul (|2QQ9f ) to analyze the Fe XIII 



202 A and Fe XIV 274 A lines. They redress the 
Alfven wave interpretation suggested bv lMcIntoshl (j2QQ9h 
and independently arrive at the same conclusion as re- 
ported in our work: the analysis of line profile asymme- 
tries shows that multiple component flows (and not only 
Alfven waves) play an importa nt role in th e plasma dy- 
namics in coronal dimmings. iPeterl (j2010h finds in an 
active region a narrow line core and a broader minor 
component blueshifted by up to 50 kms -1 for the Fe XV 
284 A line. iBrvans et al.l (|2010f ) showed that, in outflows 
observed near an active region, the Fe XII and Fe XIII 
line profiles are better represented by double- Gaussian 
fits. We find this behavior not only in upflows associated 
with dimmings events, but also in downflows associated 
with loop structures. 
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Universities of: Cambridge (UK), George Mason (USA), 
and Florence (Italy). 
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APPENDIX 

DOUBLE-GAUSSIAN FIT AND TESTS ON SYNTHETIC LINE PROFILES 

We used a double- Gaussian fit procedure including a constant and a linear term, based on a random-restart hill- 
climbing method. Basically, a classical gradient-expansion algorithm (CURVEFIT .pro in IDL) computes a non linear 
least-squares fit. It is restarted more than 10 4 times, finally selecting the solution presenting the best \ 2 - We used 
random starting estimates that could vary relatively widely: peak intensity between and the maximum of the 
spectrum in the recorded window, line center within the wavelength interval of this window and line width between 
the instrumental width and 1/4 of the spectral window width. To be selected, the solution of the gradient expansion 
had to keep into those limits. Hereafter, one has to understand the solution of a "double-Gaussian fit" as the best 
solution of multiple iterations. 

We tested the reliability of this procedure on synthetic line profiles. First we simulated a double- Gaussian profile 
with components similar to what we find in the contoured area of the dimming in data set lc: a component at rest 
and an additional one with 40% amplitude, offset by a few tens of kilometers per second, but having the same width 
as the static component. We applied the double-Gaussian fit thousands of times on the same synthetic profile (i.e., 
each solution of a fit being itself the results of thousands of iterations). We then analyzed the resulting statistics. 
Even without added statistical noise on the synthetic profiles, the amplitude, velocity and line width of the individual 
components cannot be retrieved with high precision. Not surprisingly, the parameters of the component having the 
lowest amplitude present the largest dispersion. Its amplitude is on average overestimated by 25%. The distribution 
of retrieved component centers is rather flat over an interval more than 60 kms -1 wide, with a median displaced by 
about 25 kms -1 as compared to the initial value. This means that the center of the small (and dynamic, in our case) 
component is on average artificially found redshifted as compared to the original component. We remind that an EIS 
spectral pixel is equivalent to ~ 35 kms -1 . Likewise, the line width of the small component is on average overestimated 
by 15%. In compensation, the center of the static (and more intense) component is artificially blueshifted by about 
10 kms -1 , which is similar to what we find in real data. When we add statistical noise, the results are not modified 
very much. This analysis holds for profiles where no clear double-peak appears in the resulting profile; it is of course 
easier to retrieve the correct parameters when the two components are well separated. 

In a second kind of verification, we simulated a single, large Gaussian and try to fit it with two Gaussians. Not 
surprisingly, when no noise is added, the fit provided two nearly identical Gaussians, with an amplitude half of the 
original one. When we add some statistical noise, the fitting procedure starts to "lock" on some local bumps on the 
profile and finds two components that can be separated by up to several tens of kilometers per second. On average, 
they nevertheless remain separated only by about ±10 kms -1 . 
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All these results correspond to a few particular sets of parameters and cannot be compared directly to our data set 
that more likely corresponds to a wider range of parameters. But it gives us important information about the reliability 
of results of our double- Gaussian fit: the presence of two components (at least) is real in the pixels with large non- 
thermal velocity. The fact that one of them is at rest and the other one is highly Doppler-shifted is statistically relevant, 
even though it is impossible to have confidence in the fit parameters for a given pixel. Our analysis also shows that 
one should be careful in interpreting the results of our double- Gaussian fit for pixels with small non-thermal velocity. 
Only when the coefficient of asymmetry and/or peakedness are non- negligible, one can trust the double- Gaussian fit. 

EFFECT OF REBINNING ON MEASURED DISTORTIONS 

Due to the orbital variation of the line center and to the spatial inhomogeneity of the Doppler shift, the rebinning 
(cf Sect. 14.3ft can be an additional source of distortion for the line profiles. We have to verify, then, that the rebinning 
does not contribute much to the distortion we want to emphasize. For that, we simulated a 3 x 3 rebinning but with 
the single-Gaussian fitted profiles of Fe XII instead of the real data. The derived maps of asymmetry and peakedness 
(not shown) present no particular pattern, and are absolutely not comparable to what we obtained in Fig. [3] without 
rebinning. In fact, the coefficients are nearly equal to for almost all pixels. We can then conclude that the rebinning 
does not have any effect on revealing the distortion of the line profiles qualitatively. This also demonstrates that the 
source of the distortion has nothing to do with the resolution scales: it has to be found in the LOS integration, in 
scales much smaller than the resolution scale or in temporal variation. 
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